Timely Monitoring of SARS-CoV-2 RNA Fragments in Wastewater Shows the Emergence of JN.1 (BA.2.86.1.1, Clade 23I) in Berlin, Germany.

The importance of COVID-19 surveillance from wastewater continues to grow since case-based surveillance in the general population has been scaled back world-wide. In Berlin, Germany, quantitative and genomic wastewater monitoring for SARS-CoV-2 is performed in three wastewater treatment plants (WWTP) covering 84% of the population since December 2021. The SARS-CoV-2 Omicron sublineage JN.1 (B.2.86.1.1), was first identified from wastewater on 22 October 2023 and rapidly became the dominant sublineage. This change was accompanied by a parallel and still ongoing increase in the notification-based 7-day-hospitalization incidence of COVID-19 and COVID-19 ICU utilization, indicating increasing COVID-19 activity in the (hospital-prone) population and a higher strain on the healthcare system. In retrospect, unique mutations of JN.1 could be identified in wastewater as early as September 2023 but were of unknown relevance at the time. The timely detection of new sublineages in wastewater therefore depends on the availability of new sequences from GISAID and updates to Pango lineage definitions and Nextclade. We show that genomic wastewater surveillance provides timely public health evidence on a regional level, complementing the existing indicators.

SARS-CoV-2 infection dynamics revealed by wastewater sequencing analysis and deconvolution.

The use of RNA sequencing from wastewater samples is a valuable way for estimating infection dynamics and circulating lineages of SARS-CoV-2. This approach is independent from testing individuals and can therefore become the key tool to monitor this and potentially other viruses. However, it is equally important to develop easily accessible and scalable tools which can highlight critical changes in infection rates and dynamics over time across different locations given sequencing data from wastewater. Here, we provide an analysis of lineage dynamics in Berlin and New York City using wastewater sequencing and present PiGx SARS-CoV-2, a highly reproducible computational analysis pipeline with comprehensive reports. This end-to-end pipeline includes all steps from raw data to shareable reports, additional taxonomic analysis, deconvolution and geospatial time series analyses. Using simulated datasets (in silico generated and spiked-in samples) we could demonstrate the accuracy of our pipeline calculating proportions of Variants of Concern (VOC) from environmental as well as pre-mixed samples (spiked-in). By applying our pipeline on a dataset of wastewater samples from Berlin between February 2021 and January 2022, we could reconstruct the emergence of B.1.1.7(alpha) in February/March 2021 and the replacement dynamics from B.1.617.2 (delta) to BA.1 and BA.2 (omicron) during the winter of 2021/2022. Using data from very-short-reads generated in an industrial scale setting, we could see even higher accuracy in our deconvolution. Lastly, using a targeted sequencing dataset from New York City (receptor-binding-domain (RBD) only), we could reproduce the results recovering the proportions of the so-called cryptic lineages shown in the original study. Overall our study provides an in-depth analysis reconstructing virus lineage dynamics from wastewater. While applying our tool on a wide range of different datasets (from different types of wastewater sample locations and sequenced with different methods), we show that PiGx SARS-CoV-2 can be used to identify new mutations and detect any emerging new lineages in a highly automated and scalable way. Our approach can support efforts to establish continuous monitoring and early-warning projects for detecting SARS-CoV-2 or any other pathogen.

Contrasting levels of hybridization across the two contact zones between two hedgehog species revealed by genome-wide SNP data.

Hybridization and introgression have played important roles in the history of various species, including lineage diversification and the evolution of adaptive traits. Hybridization can accelerate the development of reproductive isolation between diverging species, and thus valuable insight into the evolution of reproductive barrier formation may be gained by studying secondary contact zones. Hedgehogs of the genus Erinaceus, which are insectivores sensitive to changes in climate, are a pioneer model in Pleistocene phylogeography. The present study provides the first genome-wide SNP data regarding the Erinaceus hedgehogs species complex, offering a unique comparison of two secondary contact zones between Erinaceus europaeus and E. roumanicus. Results confirmed diversification of the genus during the Pleistocene period, and detected a new refugial lineage of E. roumanicus outside the Mediterranean basin, most likely in the Ponto-Caspian region. In the Central European zone, the level of hybridization was low, whereas in the Russian-Baltic zone, both species hybridise extensively. Asymmetrical gene flow from E. europaeus to E. roumanicus suggests that reproductive isolation varies according to the direction of the crosses in the hybrid zones. However, no loci with significantly different patterns of introgression were detected. Markedly different pre- and post-zygotic barriers, and thus diverse modes of species boundary maintenance in the two contact zones, likely exist. This pattern is probably a consequence of the different age and thus of the different stage of evolution of reproductive isolating mechanisms in each hybrid zone.

The Core Human Microbiome: Does It Exist and How Can We Find It? A Critical Review of the Concept.

The core microbiome, which refers to a set of consistent microbial features across populations, is of major interest in microbiome research and has been addressed by numerous studies. Understanding the core microbiome can help identify elements that lead to dysbiosis, and lead to treatments for microbiome-related health states. However, defining the core microbiome is a complex task at several levels. In this review, we consider the current state of core human microbiome research. We consider the knowledge that has been gained, the factors limiting our ability to achieve a reliable description of the core human microbiome, and the fields most likely to improve that ability. DNA sequencing technologies and the methods for analyzing metagenomics and amplicon data will most likely facilitate higher accuracy and resolution in describing the microbiome. However, more effort should be invested in characterizing the microbiome's interactions with its human host, including the immune system and nutrition. Other components of this holobiontic system should also be emphasized, such as fungi, protists, lower eukaryotes, viruses, and phages. Most importantly, a collaborative effort of experts in microbiology, nutrition, immunology, medicine, systems biology, bioinformatics, and machine learning is probably required to identify the traits of the core human microbiome.

Multiple types of genomic variation contribute to adaptive traits in the mustelid subfamily Guloninae.

Species of the mustelid subfamily Guloninae inhabit diverse habitats on multiple continents, and occupy a variety of ecological niches. They differ in feeding ecologies, reproductive strategies and morphological adaptations. To identify candidate loci associated with adaptations to their respective environments, we generated a de novo assembly of the tayra (Eira barbara), the earliest diverging species in the subfamily, and compared this with the genomes available for the wolverine (Gulo gulo) and the sable (Martes zibellina). Our comparative genomic analyses included searching for signs of positive selection, examining changes in gene family sizes and searching for species-specific structural variants. Among candidate loci associated with phenotypic traits, we observed many related to diet, body condition and reproduction. For example, for the tayra, which has an atypical gulonine reproductive strategy of aseasonal breeding, we observed species-specific changes in many pregnancy-related genes. For the wolverine, a circumpolar hypercarnivore that must cope with seasonal food scarcity, we observed many changes in genes associated with diet and body condition. All types of genomic variation examined (single nucleotide polymorphisms, gene family expansions, structural variants) contributed substantially to the identification of candidate loci. This argues strongly for consideration of variation other than single nucleotide polymorphisms in comparative genomics studies aiming to identify loci of adaptive significance.

The complete mitochondrial genome of the pentastomid Linguatula arctica (Pentastomida) from reindeer (Rangifer tarandus) in Northern Norway.

Here, we present the first complete mitochondrial genome of the pentastomid Linguatula arctica collected from the nasal passages of a reindeer (Rangifer tarandus) in Norway. The full length mitochondrial genome of L. arctica, which measures 14,789 bp in length, contains 13 protein-coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes. A clear A + T bias is observed in the mitogenome of L. arctica with an overall base composition of 32.6% A, 27.5% T, 32.8% C, and 7,1% G., and a GC content of 39.9%. The gene arrangement is identical to that of previously described pentastomid mitogenomes.

The developmental transcriptome for Lytechinus variegatus exhibits temporally punctuated gene expression changes.

Embryonic development is arguably the most complex process an organism undergoes during its lifetime, and understanding this complexity is best approached with a systems-level perspective. The sea urchin has become a highly valuable model organism for understanding developmental specification, morphogenesis, and evolution. As a non-chordate deuterostome, the sea urchin occupies an important evolutionary niche between protostomes and vertebrates. Lytechinus variegatus (Lv) is an Atlantic species that has been well studied, and which has provided important insights into signal transduction, patterning, and morphogenetic changes during embryonic and larval development. The Pacific species, Strongylocentrotus purpuratus (Sp), is another well-studied sea urchin, particularly for gene regulatory networks (GRNs) and cis-regulatory analyses. A well-annotated genome and transcriptome for Sp are available, but similar resources have not been developed for Lv. Here, we provide an analysis of the Lv transcriptome at 11 timepoints during embryonic and larval development. Temporal analysis suggests that the gene regulatory networks that underlie specification are well-conserved among sea urchin species. We show that the major transitions in variation of embryonic transcription divide the developmental time series into four distinct, temporally sequential phases. Our work shows that sea urchin development occurs via sequential intervals of relatively stable gene expression states that are punctuated by abrupt transitions.

Improving draft genome contiguity with reference-derived in silico mate-pair libraries.

Background: Contiguous genome assemblies are a highly valued biological resource because of the higher number of completely annotated genes and genomic elements that are usable compared to fragmented draft genomes. Nonetheless, contiguity is difficult to obtain if only low coverage data and/or only distantly related reference genome assemblies are available. Findings: In order to improve genome contiguity, we have developed Cross-Species Scaffolding-a new pipeline that imports long-range distance information directly into the de novo assembly process by constructing mate-pair libraries in silico. Conclusions: We show how genome assembly metrics and gene prediction dramatically improve with our pipeline by assembling two primate genomes solely based on ∼30x coverage of shotgun sequencing data.

The complete mitochondrial genome of the lesser sac-winged bat Saccopteryx leptura (Chiroptera: Emballonuridae) from Costa Rica.

Here we present the first complete mitochondrial genome of the lesser sac-winged bat Saccopteryx leptura (Chiroptera: Emballonuridae) from Costa Rica, assembled from next-generation sequencing data. The mitogenome of Saccopteryx leptura measures 16,577 bp in length, and contains 13 protein-coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes. A slight A + T bias was observed in the mitogenome of Saccopteryx leptura with an overall base composition of 31.5% A, 28.3% T, 25.8% C, and 14.2% G, and a GC content of 40.1%. The gene arrangement was identical to that of previously described bat mitogenomes.

The complete mitochondrial genome of the pentastomid Armillifer grandis (Pentastomida) from the Democratic Republic of Congo.

We present the first complete mitochondrial genome of the pentastomid Armillifer grandis (Arthropoda: Pentastomida) collected from the lungs of a rhinoceros viper (Bitis nasicornis) in the Democratic Republic of Congo. The full length mitochondrial genome of Armillifer grandis, which measures 16,073 bp in length, contains 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. A clear A + T bias is observed in the mitogenome of Armillifer grandis with an overall base composition of 34.6% A, 29.4% T, 29% C, and 6.9% G, and a GC content of 35.9%. The gene arrangement is identical to that of previously described pentastomid mitogenomes.

The complete mitochondrial transcript of the red tail loach Yasuhikotakia modesta as assembled from RNAseq (Teleostei: Botiidae).

Here, we present the first complete mitochondrial genome of the Red Tailed Loach Yasuhikotakia modesta (Teleostei: Botiidae) from Thailand, assembled from next-generation transcriptome sequencing data. The assembled transcript corresponds to the full length mitochondrial genome of Y. modesta, which measured 16,865 bp in length, and contained 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. A slight A + T bias was observed in the mitogenome of Y. modesta with an overall base composition of 32.2% A, 25.8% T, 26.4% C, and 15.4% G, and a GC content of 41.8%. The gene arrangement was identical to that of previously described loach mitogenomes.

The complete mitochondrial genome of the viviparous freshwater snail Tylomelania sarasinorum (Caenogastropoda: Cerithioidea).

Here, we present the first complete mitochondrial genome within the gastropod family Pachychilidae, using the viviparous freshwater snail Tylomelania sarasinorum. This species is a representative member of the lacustrine Tylomelania radiations of the Malili-Lakes-System (Sulawesi, Indonesia). The mitochondrial genome was 16,632 bp long and contained 13 protein-coding genes, 2 rRNA genes and 22 tRNA genes. A pronounced A + T bias was observed with an overall base composition of 29.5% A, 35.7% T, 18.3% G and 16.6% C. Tylomelania sarasinorum exhibited a novel mitochondrial gene arrangement, differing from all Caenogastropoda mitochondrial genomes published to date.

LTR retroelements are intrinsic components of transcriptional networks in frogs.

BACKGROUND: LTR retroelements (LTR REs) constitute a major group of transposable elements widely distributed in eukaryotic genomes. Through their own mechanism of retrotranscription LTR REs enrich the genomic landscape by providing genetic variability, thus contributing to genome structure and organization. Nonetheless, transcriptomic activity of LTR REs still remains an obscure domain within cell, developmental, and organism biology. RESULTS: Here we present a first comparative analysis of LTR REs for anuran amphibians based on a full depth coverage transcriptome of the European pool frog, Pelophylax lessonae, the genome of the African clawed frog, Silurana tropicalis (release v7.1), and additional transcriptomes of S. tropicalis and Cyclorana alboguttata. We identified over 1000 copies of LTR REs from all four families (Bel/Pao, Ty1/Copia, Ty3/Gypsy, Retroviridae) in the genome of S. tropicalis and discovered transcripts of several of these elements in all RNA-seq datasets analyzed. Elements of the Ty3/Gypsy family were most active, especially Amn-san elements, which accounted for approximately 0.27% of the genome in Silurana. Some elements exhibited tissue specific expression patterns, for example Hydra1.1 and MuERV-like elements in Pelophylax. In S. tropicalis considerable transcription of LTR REs was observed during embryogenesis as soon as the embryonic genome became activated, i.e. at midblastula transition. In the course of embryonic development the spectrum of transcribed LTR REs changed; during gastrulation and neurulation MuERV-like and SnRV like retroviruses were abundantly transcribed while during organogenesis transcripts of the XEN1 retroviruses became much more active. CONCLUSIONS: The differential expression of LTR REs during embryogenesis in concert with their tissue-specificity and the protein domains they encode are evidence for the functional roles these elements play as integrative parts of complex regulatory networks. Our results support the meanwhile widely accepted concept that retroelements are not simple "junk DNA" or "harmful genomic parasites" but essential components of the transcriptomic machinery in vertebrates.